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of electromagnetic backgrounds based on the worldline formalism. As 
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Introduction 



The worldline formalism was invented by Feynman simultaneously with 
modern relativistic second-quantized QED, but for a long time it was used 
only occasionally for actual calculations. Eventually, the observation of a 
close relation between the worldline formalism and the infinite-string-tension 
limit of string path integrals triggered further developments of the world- 
line formalism for QCD || and QED ||, particularly for gauge particle 
amplitudes. Certain computational advantages of this formalism were sub- 
sequently recognized and led to numerous applications. Among them, the 
progress achieved for QED amplitudes with all-order couplings to an external 
background field is particularly remarkable ||, |5], ||. The formalism works 
most elegantly for constant electromagnetic background fields and can be 
extended to a derivative expansion in the electromagnetic background |7|, § . 
For a comprehensive review of the worldline formalism in QED and beyond, 
see 0. 

In the present work, we intend to demonstrate that the worldline formal- 
ism is moreover perfectly suited for numerical computations. This is because 
the path integral over closed loops in spacetime can be approximated by a 
finite ensemble of loops, which allows for a simple and fast evaluation of ex- 
pectation values. The latter include observables depending on an arbitrary 
background field. 

We illustrate this proposal by means of a simple example: the one-loop 
effective action of Euclidean scalar QED in three dimensions. The compu- 
tational task boils down to the calculation of the Wilson loop expectation 
value using a loop ensemble. 

As a first step, we verify the method in Sect. |2]by considering the constant- 
magnetic-field case, which can also be solved analytically. This serves as a 
benchmark test for the concrete procedure that we propose. In a second 
step, the numerical method is applied to a magnetic field resembling a step 
function in one spatial direction (Sect. |3|). This idealized field configuration 
represents the simplest configuration, which is inaccessible to the standard 
analytical method for inhomogeneous fields: the derivative expansion. As 
our main result, we observe diffusion of the magnetic field: the field takes 
influence on the region of vanishing background field by inducing a nonzero 
action density therein. 

Our conclusions are summarized in Sect. [|, where possible generalization 
and the road to further application of worldline numerics are sketched. 
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1 The worldline approach 
to functional determinants 



1.1 Setup 

Our starting point is the unrenormalized Euclidean one-loop effective action 
of scalar QED in D dimensions in worldline representation || , corresponding 
to the determinant of the gauge-covariant Klein-Gordon operator 

T [A] = / -^-e~ AT / Vx(r)e o ^ 4 "\ (1) 

l/A 2 x(T)=x(0) 

where the superscript "1" indicates the one-loop levelQ; a gauge-invariant UV 
regularization at a scale A has been performed at the lower bound of the T 
integral for the sake of defmiteness. In Eq. ([J), we encounter a path integral 
over closed loops in spacetime. Note that there are no other constraints 
to the loops except differentiability and closeness; in particular, they can be 
arbitrarily self-intersecting and knotty. The normalization can be determined 
from the zero- field limit, 

Solving Eq. (Q) for M and inserting this into Eq. (JT|) leads us to the compact 
formula 

oo 

r V] = (4^5572 / dD *« I r^ e ' m2T{W[A]) - (3) 

1/A 2 

Here we have split off the integral over the zero-modes of the path integral, 
f d D XQ, where xq, the so-called loop center of mass, corresponds to the aver- 
age position of the loop: Xq := (1/T) J Q drx^ij). In Eq. (^|), we introduced 
the Wilson loop 

T 

— ie f dr x(t)-A(x(t)) . „ , ... 

W[A(x)\ = e = e- ie $ dx - A(x \ (4) 

1 In an abuse of nomenclature, we call the one-loop contribution to the total effective 
action also "effective action" in the following. The reader should always keep in mind that 
the Maxwell term (and all higher-order terms) must be added. 
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and ((. . . )) x denotes the expectation value of ( . . . ) evaluated over an en- 
semble of x loops; the loops are centered upon a common average position 
Xq ( "center of mass" ) and are distributed according to the Gaussian weight 
exp[-/ T drf ]. 

The following r-integral substitution, r =: Tt, is of crucial importance 
for the numerical realization of the path integral; it suggests introducing unit 
loops y, 

y(t):=^=x(Tt), te[0,l], (5) 



which are parameterized with a unit propertime t. The remaining integrals 
can be rewritten accordingly: 

T 2 1 2 T 1 

dr = [alt ^p, / oItx{t) ■ A{x{t)) = [ dty(t) ■ [VTA(VTy(t))}. 
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(6) 



The important advantage is constituted by the fact that the expectation 
value of W[A] can now be evaluated over the unit-loop ensemble y, 

(W[A(x)}) x = (W[VTA(VTy)]) y , (7) 

where the exterior T-propertime dependence occurs only as a scaling factor 
of the gauge field and its argument. In other words, while approximating the 
loop path integral by a finite ensemble of loops, it suffices to have one single 
unit-loop ensemble at our disposal; we do not have to generate a new loop 
ensemble whenever we go over to a new value of T. 

Inserting Eq. (|7|) into Eq. (|j), we arrive at the final formula (dropping 
the subscript of x ): 

oo 

r'l-4] = ^/A / ^-^{w^A^y) 

1/A 2 ^ 



at., 



where we have formally added counter-terms (c.t.) which correspond to a 
renormalization of the physical parameters. The details of the c.t.'s depend 
on the dimension, but not on the type of the background field, and will be 
discussed in the following. 
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1.2 Renormalization 



For an accurate evaluation of the effective action of Eq. (|j), an analytic 
calculation of the counter-terms is inevitable in order to obtain accurate 
results by the numerical procedure. In the following, we will confine ourselves 
to the important cases of D = 3 and D = 4 dimensions at the one-loop level. 

Let us first discuss the analytical treatment: while for D = 3 the effec- 
tive action is rendered finite by dropping a field-independent constant, the 
complete effective action in the 4-dimensional case is renormalized in such a 
way that the quadratic term in the field strength is given by 



1 



quadr. 



d x F^ v (x) F^x) 



(9) 



where F fJLy (x) is the field strength tensor, and represents the renormalized 
coupling. In practice, this is done by subtracting the (infinite) quadratic 
term of the one-loop contribution T 1 and absorbing it into the bare Maxwell 
term. 

These UV divergencies of the effective action emerge from the lower bound 
of the propertime integral in Eq. i.e., T — > 0. For small values of the 
propertime, the Wilson loop expectation value can be calculated exactly (e.g., 
using heat-kernel techniques): 



W[VTA(Vfy)} 



1 " -T 2 F, v [A]{x) F, u [A]{x) + 0(T 4 



(10) 



where x is the loop center of mass. The key observation is that, even for 
D = 4, the terms of order T 4 are UV-finite upon the propertime integration 
in Eq. (|8|). Hence, the completely renormalized effective action, which is 
suitable for numerical simulations, is 



r^A] 



i 



where 



c 3 



(4tt) c / 2 



1 1 

967r m 



d D x 



dT 



-m 2 T 



T(D/2)+l 



+c D J d D x [A] (x) F^ [A] (x) 



+ 0(1/ A), c 4 



W[VfA(Vfy)] - l) (11 
/ y 

+ ^T 2 F, I/ [A}(x)F, u [A}(x 



12 16tt 2 



A 2 

In— -C 



0(m 2 /A 2 



(12) 
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Here A denotes the UV cutoff, and C is Euler's number. In the D = 4 case, 
the term ~ C4 will be absorbed into the bare Maxwell term, defining the 
running of the coupling, and the cutoff can subsequently be sent to infinity. 
In D = 3, the theory is super- renormalizable and the term ~ C3 represents 
only a finite shift of the coupling, introducing no running. In the first two 
lines of Eq. (|TT|) , we have already sent the cutoff A two infinity, leaving us with 
a perfectly finite expression. Note that in the case D = 3, the term of order 
T 2 in Eq. ( |T0D corresponds to a singularity 1/y/T, which is integrable. Hence, 
the subtraction of this term from the propertime integral is not mandatory. 
However, aiming at worldline numerics, the perfect control over the behavior 
of the integrand at small propertime T is required in order to augment the 
accuracy of the numerical evaluation of the propertime integral. 

Now, the numerical renormalization is similar in spirit to the analytical 
one, but is complicated by a further problem: evaluating (W) with the aid of 
the loop ensemble does not produce the small-T behavior of Eq. (|T0D exactly, 
but, of course, only within the numerical accuracy. Unfortunately, even 
the smallest deviation from Eq. (|T0D will lead to huge errors, if we naively 
plug such a result into Eq. (|12"D ; this is because of the factors of T in the 
denominator of the integrand for T — > 0. 

Our solution to this problem is to fit the numerical result for (W) to a 
polynomial in T in the vicinity of T = 0, employing Eq. fllQ|) as a constraint 
for the first coefficients. Of course, such a fit is completely justified because 
of our exact knowledge of (W) for small T. This fit not only represents the 
renormalization procedure, solving the UV problems, but also facilitates a 
more precise estimate of the error bars (see below). Finally, employing this 
fitting procedure only close to T = 0, the infrared behavior (T — > 00) of the 
integrand remains untouched, and our approach is immediately applicable, 
also in the case m — 0. 

Both renormalization procedures, the analytical as well as the numerical, 
generalize straightforwardly to higher dimensions; here, either additional sub- 
tractions from the integrand (for the analytical case), or polynomial fits with 
higher-order constraints (for the numerical case) are required: for example, 
increasing the dimension by 2 requires one more subtraction/constraint in 
the small-T behavior of the integrand. 

Finally, it should not be concealed that these procedures become increas- 
ingly difficult to higher order in perturbation theory for massive theories in 
D > 4. Then, a mass renormalization is also necessary, requiring careful 



analyses of the UV behavior of double propertime integrals [12 
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1.3 Numerical simulation 



Now, the route to the effective action is clear: 

(1) generate a unit loop ensemble distributed according to the weight 
exp[— J" 1 dty 2 /A\, e.g., employing the technique of normal (Gaussian) 
deviates; 

(2) compute the integrand for arbitrary values of T (and x); this involves 
the evaluation of the Wilson loop expectation value for a given back- 
ground gauge field; 

(3) perform the renormalization procedure and add the c.t.'s to the inte- 
grand; 

(4) integrate over the propertime T in order to obtain the Lagrangian, and 
also over x for the action. 

From a general viewpoint, there are two sources of systematical error 
which are introduced by reducing the degrees of freedom from an infinite 
to a finite amount: first, the loop path integral has to be approximated by 
a finite number of loops; second, the propertime t of each loop has to be 
discretized. Contrary to this, the spacetime does not require discretization, 
i.e., the loop ensemble is generated in the continuum. 

Of course, the various steps can be carried out using different numerical 
methods; let us outline our choice of tools in detail. For this purpose, let 
Ni denote the number of loops which are used to estimate the Wilson loop 
expectation value in Eq. (^), and ni the number of space points which are 
employed to specify a particular loop. The points of a particular loop y^ % = 
1 . . . rii, are generated by a standard heat bath algorithm, where the boundary 
conditions y\ = 0, y n i = are enforced. After a proper thermalization, all 
coordinates y^ are shifted equally in order to implement the center of mass 
condition (y 4 ) = x . This procedure is then repeated Ni times to generate 
the loop ensemble. 

Approximating (W) of Eq. (|7|) by an average over a finite number of 
loops, the standard deviation provides an estimate of the statistical error. 
Approximating the loops by the finite number ni of space points results in 
a systematic error that can be estimated by repeating the calculation for 
several values n\. The number rii should be chosen large enough to reduce 
this systematic error to well below the statistical one. It will turn out that 
the choice Nl = 1000 and rii = 100 yields results, for the applications below, 
which are accurate at the per cent level. We stress, however, that ~ 20000 
(dummy) heat bath steps are required to properly thermalize the loop en- 
semble. 
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2 Benchmark test: 

constant magnetic background field 

Since analytic results are available for the case of a constant magnetic back- 
ground field B, we will investigate in this section the efficiency of our nu- 
merical loop approach to the scalar functional determinant .Q For simplicity, 
we consider three Euclidean spacetime dimensions, D = 3. Up to a field- 
independent constant, the one-loop effective-action density (Lagrangian) C 1 
for scalar QED in D = 3 is given by M, ITT 



C 1 = Y l {B)/V z 



£? 3 / 2 /m 2 
(4tt) 3 / 2 9 \~B 



(13) 



where V* denotes the volume and 



9(z) 



dT 



I(T, z) 



I(T, z) 



1 

J^2 



T 



sinhT 



- 1 



exp{— zT} 



(14) 



The exact integrand J(T, z) in Eq. (|IJ|) is compared with our loop estimate 
in Fig. [I] for the case m 2 = 0, i.e., z = 0. The numerical estimate for g(m 2 /B) 
is also shown. The agreement between the two curves is satisfactory and the 
exact results for I(T, z) as well as for g(m 2 /B) lie well within the error bars 
produced by the numerics. 

As mentioned above, the error bars correspond to the statistical error of 
the ensemble average. Regarding the quantity (W) with its dependence on 
T, the original error bars are rather independent of T; but multiplying (W) 
by a T-dependent function (cf. Eq. (H)) causes a modulation of the error 
bars. In particular for T — > 0, the 1/y/T singularity leads to an unbounded 
enhancement of the error bars for small T in the function I(T,x) (see left 
panel of Fig. [I]). Fortunately, the integrand in this regime is known exactly as 
given in Eq. (|irj) in the form of an asymptotic series. As mentioned above for 
the numerical renormalization procedure, this information can be used for 
a constraint fit of the numerical result to a polynomial in T around T = 0, 
keeping the first terms corresponding to Eq. ( ]TTJ| ) fixed. Then, the error bars 
of the ensemble average translate into errors for the higher coefficients of the 
polynomial. This polynomial is then matched to the pure numerical result 
at that value of T where the error bars of the two results are comparable. 
The final result is visible on the left panel of Fig. [I] 

incidentally, another exactly solvable system is known with a background field 
resembling a solitonic profile in one dimension. 



s 




T m 2 /B 



Figure 1: Propertime integrand I(T,z = 0) (left panel) and integral g(m 2 /B) 
(right panel) of the one-loop effective action for the case of a constant magnetic 
background field. The analytically known exact results (solid lines) are compared 
with the numerical findings (circles with error bars). 

It should be noted, that the error bars in Fig. [1] (and the following figures) 
are highly correlated from point to point, since the same loop ensemble has 
been used for the evaluation of each point. This correlation can, of course, 
easily be reduced by updating the loop ensemble with a few heat bath steps 
inbetween at the expense of computational time. 

3 Magnetic-field diffusion 

In this section, we study the one-loop effective action of Eq. ( |TT| ) for the 
case of a magnetic background field, resembling a step function in space. 
In particular, we consider a time-like constant background field B, i.e., a 
field which is independent of the third coordinate called Euclidean time. We 
choose the B field, being a (pseudo-) scalar over the spatial xy plane, as 

B{x, y) = -9(x) B , A{x, y) = 6{x) X - {y, -x) B , (15) 

where 9(x) is the step function, i.e., 

9(x) = 1 , Va; > , 6(x) = , Vx < . 

Note that because of the sudden variation of the background field at the 
step, the effective action cannot be obtained within a derivative expansion. 
By contrast, such a discontinuity represents no obstacle for the numerical 
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worldline approach proposed in this work. Discontinuities do not induce 
(artificial) singularities, but are smoothly controlled by the properties of the 
loop ensemble. This ensemble resembles a cloud, exhibiting finite extension 
and slowly varying density, and being centered at some particular spacetime 
point xq which is the center of mass of each loop in the cloud. While running 
with this xo towards and across the step, that part of the volume of the 
loop cloud which "feels" the magnetic field increases smoothly; therefore, 
the effective-action density (effective Lagrangian) C 1 , being the propertime 
integrated information of what the loop cloud measures, will be smooth, too. 

In order to numerically evaluate the effective action for the case of Eq. ([TBI), 
we employ the loop approach outlined in the previous section. This approach 
provides for statistical error bars which allow to estimate whether the number 
of loop ensembles is large enough. In order to study the systematic errors, 
such as the limited number of spacetime points ni specifying a loop and the 
limited number of thermalization sweeps, we shall employ three sets of 
loop ensembles (see Table 1). Defining 

g xB 1 ^ := C\x = (x, 0, 0)) (16) 



in analogy to Eq. ( J13|) , the final numerical result for the effective action 
is shown in Fig. |2|, left panel. As expected, the effective-action density is 
nonzero even in the region x < where the background field B(x) vanishes. 
For increasing mass, the contributions from large propertime T are exponen- 
tially suppressed in the integrand. Since the propertime controls the size of 
the loop cloud, the large loop clouds contribute less to the action density 
when the mass is large. Hence, the effective-action density for an increasing 
mass becomes reduced. (From an alternative viewpoint, the limit of large 
mass and small field are identical for dimensional reasons.) 

We find that the effective action decreases exponentially in the "forbid- 
den" region, which is depicted in Fig. 0, right panel; here we defined the 



Table 1: Simulation parameters 





ni 


Ht 


set A 


100 


150000 


set B 


75 


50000 


set C 


50 


50000 



10 



CQ 

x 

o 

CQ 

OJ 



-0.2 - 



-0.4 - 



-0.6 - 



-0.1 



-2 




o m=0, set A 
a m=0, set B 

m=0, set C 
a m=0.5 B 1 ' 2 , 
< m=0.5 B n 1 ' 2 , 

m=0.5 B„ 1 ' 2 , 



set A 
set B 
setC 



o m=0, set C 
□ m 2 =B , set C 
o m 2 =2B , set C 
a m=0, set B 
< m 2 =B , set B 
v m s =2B„, set B 
> m=0, set A 
+ m 2 =B , set A 
- m s =2B„, set A 



-Kill 



x B n 



xB„ 



Figure 2: Effective-action density in the vicinity of the magnetic step {x = 0). 
On the left panel, the diffusion profile is visualized. The logarithmic plot on the 
right panel reveals the exponential nature of the diffusion depth. 



quantity 



-// ( '-^-,xBl /2 ) =: exp -S (xBl 12 



(17) 



and plotted S(xB^ 2 ) for several values of the scalar mass m. This phe- 
nomenon of obtaining a nonzero effective-action density even in a region 
where the background field vanishes may be called quantum diffusion of the 
magnetic field. 

Large values of S imply that the effective-action density ~ g (^m 2 /B , xB^ 2 

is small due to cancellations. In this sufficiently large number of loop 

points rii and of thermalization sweeps is requested. The results indicate that 

1 /2 i 

the function S(xB ' ) can be fitted for Xy/Bo by the ansatz 

S(y) = a + /3y, y = xB x J 2 \ (18) 

where the quantities a and (3 depend on m 2 /B . The diffusion depth I of 
the magnetic field can be defined as the inverse of /3, i.e., I = \/{(3\/B~q). An 
estimate of this parameter has been plotted versus the mass m in Fig. |. As 
expected, the diffusion coefficient /3 rises with increasing mass scale m. The 
function 

2\ 1/2 



(3(m 2 /B ) = 0.7627 + 3.255 ( — ) (19) 
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Figure 3: Inverse diffusion depth (3 versus mass-to-magnetic-field ratio. The plot 
symbols depict the numerical results, whereas the solid line represents the fit of 
Eq. (19) adjusted to the optimized loop ensemble A. 



nicely fits the high quality numerical data set A (see figure §). 

In order to get an understanding of these numbers, let us perform the 
following heuristic consideration: at a first glance, one might expect that the 
magnetic diffusion obeys the simple law ~ e -const.-ma^ s i nce ^ ne mass seems 
to be the only scale in the field-free region of space; this would correspond to 
(3 = const. • (m 2 / Bq) ' 5 ■ However, if this simple law were true, the massless 
limit m = would be obscure, since then the magnetic field would diffuse 
into the field-free region without any damping. Hence, there must be an 
additional dependence of the exponent on the magnetic field in order to 
account for a reasonable massless limit. Moreover, the mass is indeed not 
the only scale in the field-free region, because the loop cloud is a nonlocal 
object that always "feels" the strength of the magnetic field. In fact, it is the 
constant first term in f3 in Eq. ( |T§D that exactly accounts for this dependence 
of the magnetic diffusion on the strength of the magnetic field. Only for larger 
masses (or weaker fields), m 2 ~ B /4, the intuitively expected diffusion law 
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of the form ~ e const - mx ' begins to dominate; in this regime, we find 

/ « 0.31/m. (20) 

In the present work, the precision of the numerics for the action density in 
the field- free region restricts the investigation to mass values of m 2 < 1.5B ; 
beyond this, the strong exponential decrease beyond the step prohibits a 
reliable analysis of the diffusion depth. 



4 Conclusion and Outlook 

Beyond any particular result of the present work, we would like to remark 
in the first place that our approach to the worldline formula (see Eq. (§)) 
for the one-loop effective action offers a vivid picture of the quantum world. 
Consider a spacetime point x at a propertime T; then, the loop ensemble 
is centered upon this point x and resembles a loop cloud with Gaussian 
"density" and "spread" . Increasing or lowering the propertime T corresponds 
to bloating or scaling down the loop cloud or, alternatively, zooming out of or 
into the microscopic world. The effective-action density at each point x finally 
receives contributions from every point of the loop cloud according to its 
Gaussian weight (times the mass term and other factors) and averaged over 
the propertime. This gives rise to the inherent nonlocality and nonlinearity 
of the effective action, because every point x is influenced by the field of any 
other point in spacetime experienced by the loop cloud. 

In particular, we considered the one-loop contribution to the effective 
action in scalar electrodynamics in three spacetime dimensions; in the be- 
ginning, we were able to reproduce the analytically well-known case of a 
constant magnetic background field, serving as a testing ground for our nu- 
merical procedures including renormalization to one-loop order. Incidentally, 
the zero-mass limit (or, alternatively, the ultra-strong magnetic field limit) 
is also covered by our approach without additional difficulties. 

We furthermore tackled the problem of a step-function-like magnetic 
background field, illustrating the stability of our approach also for discontinu- 
ous field configurations. For this case, we observed a diffusion of the magnetic 
field, i.e., nonzero action density even at a distance from the magnetic field. 
This phenomenon is obviously nonlocal, since an expansion around a point 
of zero background field gives a zero result to any finite order. But the dif- 
fusion phenomenon appears to be also nonperturbative in the same sense as 
the Schwinger mechanism of pair production |TTfl, at least for not too large 



values of the mass; this is suggested by the functional form of the diffusion 
for small mass: exp{— S(xB ' )) ~ exp(— 0.7627 • y/Box); this result cannot 
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be expanded in terms of the coupling constant, being rescaled in the field 
(only an expansion in terms of the square root of the coupling constant is 
possible). 

Perhaps the main advantage of the numerical worldline approach to func- 
tional determinants is that no a-priori information, such as certain symme- 
tries of the background field or a suitably chosen set of base functions, has 
to be exploited; if the loop ensemble has been thermalized properly, any 
background field can immediately be plugged into the algorithm. 

One drawback of the numerical approach lies in the fact that it applies 
only to Euclidean quantum field theory; this is because the action governing 
the distribution of the loops must be positive. 

The present paper paves the way to further generalizations such as fer- 
mionic functional determinants; our approach could facilitate a systematic 
numerical investigation of these determinants. Results could be compared 
with several results known from analytical considerations (see, e.g., [p~3|j ). 
Although the treatment of fermions in the worldline formalism can be most 
elegantly formulated via Grassmann representations (see, e.g., [1^|), a numer- 
ical approach has to rely on a bosonic representation of the path integrals 
PI [15]]. In practice, this means that the Dirac algebraic elements in the 
worldline action are accompanied by a path-ordering prescription. Similar 
complications occur for nonabelian gauge fields and color-charged fermions 
or scalars. Whereas such path ordering is difficult to deal with in analyt- 
ical approaches, a numerical evaluation can immediately take care of such 
a prescription. This is because the loops are discretized in the propertime 
parameter anyway, consisting of ni "links" ; the path ordering then is nothing 
but a simple (matrix-) multiplication of these links. 

Recently fl6l , a new approach was designed to improve the still un- 
pleasant situation [I7j when lattice QCD is studied at finite baryon densities: 



the basic idea to circumvent the problems |L8| with lattice fermions is a cal- 
culation of the continuum fermion determinant for arbitrary entries of the 
gluon field, which subsequently is the subject of a lattice discretization. This 
program was successfully applied to the case of heavy quarks |TB[ . The for- 
mulation of the present paper stirs the hope that the concept of [|16[ might 
be extended to the realistic case of light quarks. 

Further possible and straightforward applications can be found in the 
context of estimating quantum energies to solitonic fields (see [0), or in 
the context of thermal field theory or Casimir vacua (for analytical worldline 
approaches, see, e.g., [0]). The latter cases are connected with a compacti- 
fication of the Euclidean spacetime manifold; these topological modifications 
can be imposed directly on the loop ensemble in our approach. For example, 
compactifying the Euclidean time direction at finite temperature will result 



14 



in closed loops that wind several times around the time-like cylinder. Work 
in these directions is in progress. 

Finally, it is obvious that the step-like magnetic field is not physical at 
all in a strict sense, because such a sharp drop-off cannot be produced by 
real laboratory magnets; nevertheless, the step-like field can be regarded as 
a limiting case of a real physical situation, being useful for an estimate of 
possible effects caused by rapid variations of a magnetic field. For instance, as 
a matter of principle, magnetic diffusion can be included in the discussion of 
the Aharonov-Bohm effect: if the propagating particle is no longer considered 
as a quantum mechanical point particle, but as a quantum field, it will be 
affected by the quantum diffusion of the solenoid's magnetic field; but a 
measurable effect can only arise, if the particle's distance from the solenoid 
is of the order of a few Compton wavelengths. 

Beyond that, we would like to mention that our results for the step-like 
magnetic field are of some significance for the experiments being currently 
performed at PVLAS (Legnaro) and BMV (Toulouse) |j21 |, in which is mea- 
sured the optical birefringence of the quantum vacuum exposed to a magnetic 
field. For example in the PVLAS experiment, the polarization axis of a laser 
beam is affected by a (6-9 T)-magnetic field with a diameter of 1 m; the 
magnetic field drops off over a distance of 10 cm, and the question arises as 
to whether this drop-off will influence the rotation of the polarization axis as 
predicted by a constant-field QED calculation. 

The refractive indices of the modified vacuum are proportional to the 
energy density induced by the magnetic field (for a review, see ||22|| ). For 
weak magnetic fields, the energy density is directly related to the (Euclidean) 
effective-action density. Provided that our present results also hold for D = 4 
spinor QED at least qualitatively, we can exclude any further influence of the 
drop-off region, since the magnetic diffusion in this case (m ^> B) occurs at 
the order of a Compton wavelength 1/m. This "zero-result" is supported by 
considerations within a derivative expansion, where the natural expansion 
parameter is given by the Compton wavelength over the length of the field 
variation (e.g., ~4x 10 -12 for the PVLAS). Only the constant-field result 
integrated over the size of the magnet including the drop-off region need be 
taken into account. 
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